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ABSTRACT 

Context. The link between the duration of GRBs and the nature of their progenitors remains disputed. Short bursts (with durations of 
less than ~2 s) are less frequently observed, technically more difficult to localize, and exhibit significantly fainter afterglows. 
Aims. It is of critical importance to establish whether the burst duration can reliably distinguish the different GRB population models 
of collapsars and compact stellar mergers. The Swift GRB 090426 provides an unique opportunity to address this question. Its duration 
(Tgo = 1.28 s) places GRB 090426 firmly in the short burst population, while the high redshift (z = 2.609), host galaxy properties, 
and prompt emission spectral characteristics are more similar to those of long-duration GRBs. 

Methods. On the basis of data obtained with the Tautenburg 2m telescope (Germany) and the 7-channel imager GROND (La Silla, 
Chile), we compiled the most finely sampled light curve available for a short burst optical/NIR afterglow. The light curve was then 
analysed in a standard fashion. GROND and XRT data were used to determine the broad-band spectral energy distribution of the 
afterglow across more than three orders of magnitude. 

Results. Our data show that a light curve break exists at 0.4 days, which is followed by a steep decay. This light curve decay is 
achromatic in the optical/NIR bands, and interpreted as a post-jet break phase. The X-ray data do not disagree with this interpretation. 
Conclusions. The half-opening angle of the suspected jet as well as the luminosity of the optical afterglow provide additional evidence 
that GRB 090426 is probably linked to the death of a massive star rather than to the merger of two compact objects. 

Key words. Gamma rays: bursts - ISM: jets and outflows: individual: GRB 090426 
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C 1 . Introduction optical detections and about one third h ave redshift d etermina- 

tions based on host galaxy spectroscopy (Berger 2010). 
It is commonly accepted that long GRBs are linked to the GRB 090426 ifj an outstanding short burst (Tc)0 „ L28 s)> 
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core-collapse of massive stars (so-called collapsar events; beC ause it has by far the highest redshift known among the short 

I Woosley&Blooml |20Q6j) residing in star-forming galaxies, burst Je {z = 2 6Q9; Antonelli et aL 2009; Levesaue et al. 

while short bursts are linked to compact stellar mergers in all |pTJ AU othe r short b ursts with secure redshift measurements 

H ■ morphological types of galaxies ( |Fongetal , 2010; Nakar2007). haye z < u (Berger 2010) The redshift of GRB Q9Q426 ifj 

Since the launch of the Swift satellite ^Gehrels et al. | |2004|), therefore in much closer agree ment with the distribution of long 

about 100 long GRBs have been rapidly localized per year GRB redshifts than wkh ^ of the short burst sample Seyeral 

( |Gehrelsetal. || 2009|) . Nearly half of them have a detected op- a nts have been t forward that GRB 090426 is not re _ 

tical afterglow and one third have redshift determinations (see Med tQ ^ CQ objectfj Antonell y^ ^ 09) found 

J. Greiner s web-pagdj. Compared to the long burst sample, that ±e QRB spectral and 

ener gy properties are more compara- 

short bursts are less frequently observed, and generally followed We tQ those of coll events . LevesauHFd]<|2010) noted that 

by on average significantly fainter and less luminous afterglows the Wue star _ forming host of this burst might also be consiste nt 

(Ka nn et alJ | 2008t | 2010| ; |Nysewander et al|2009|) . By the end of wkh a collapsar origin . similarly, Xin et al. (2010) argued that 

2010, about three dozen short GRBs had been localized by Swift the deduced lower Hmit tQ the circumburst ticle number den . 

(fewer than 10 events per year). Among them, about 50% have sky (about 1Q cm - 3) is much higher than expected for a merging 

stellar system, being more characteristic of a star-forming re- 

* Present address: American River College, Department of Physics 8 lon - 

and Astronomy, 4700 College Oak Drive, Sacramento, CA 95841, USA We present additional multi-color photometry of the opti- 

1 http://www.mpe.mpg.de/~jcg/grbgen.html cal/NIR afterglow of GRB 090426 from about 0.3 to 2.5 days 
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Fig. 1. Finding chart of the afterglow of GRB 090426. The 
r'-band GROND image is a combination of OBs 6 to 10 
(Table |A~TT ). The field of view is approximately 2' x 3'. 



after the burst, showing that the optical light curve has a well- 
defined break at late times0 



2. Observations and data reduction 

Observations of the optical/NIR afterglow of GRB 090426 
were performed with the 2m telescope of the Thtiringer 
Landessternwarte Tautenburg ( TLS, Germany) and the 7-band 
multichannel imager GROND dGreiner et al . 2008) mounted at 
the 2.2m ESO/MPI telescope on La Silla (Chile). 

Tautenburg started observing in the I c and R c band at 19:58 
UT on 26 April 2009 and stayed on the field for 90 min. The av- 
erage airmass was 1.1 and the average seeing 1.2 arcsec. About 
3.5 hr after the end of the TLS observations, at 01:08 UT on 
27 April 2009 (12.3 hr after the GRB trigger), GROND started 
following the afterglow once the target became visible over La 
Silla. Observations continued until 04:55 UT at an average see- 
ing of 1.2 arcsec and an average airmass of 2.5, during which 
GRO ND w as able to detect the afterglow in 10 different Ob£| 
(Table IaTT i. Second and third epoch observations were obtained 
the following two nights. 

Data reduction was performed in a standard fashion. TLS 
and GROND data were analysed through standard PSF pho- 
tometr y usi ng DAOPHOT and ALLSTAR tasks of IRAF 
(lTodvlll9 93). in a s imil ar way to the proc edure described in 
iKriihler et all (120081) and lYoldas et al l (120081) . Calibrations were 
performed against the SDSS, using the transformation equa- 
tions of LuptorQ for the TLS data. Magnitudes were corrected 
for Galactic extincti on, assuming E(B — V) = 0.017 mag 
dSchlegel et al.|[l998l) and a ratio of total-to-selective extinction 
of fly = 3.1. 



3. Results 

3.1. Spectral energy distribution (SED) 

The host galaxy of GRB 090426 is an extended source (about 
2") and the afterglow lies above the brightest part of this 



800 1000 
Wavelength [nm] 

Fig. 2. The SED of the afterglow at a mean time of f=53 ks 
(Table lA2b . The data is corrected for Galactic extinction and un- 
derlying host galaxy contribution. The fit is a single power law 
(X 2 /d.o.f. = 0.17). There is no evidence of extinction by dust in 
the GRB host galaxy (Ay(host)=0 mag). For the fit, the redshift 
was fixed to z = 2.61, the Lyman a absorption affects the g' band 
slightly. From left to right, we present results for the g'r'i'z' JHK 
bands. 



galaxy, the N-E knot (lAntonelli et al J 120091) . For the galaxy as 
well as this knot, g'r'i 'z' magnitudes were previously published 
(lAntonelli et al . 2009), but in the NIR bands only upper limits 
are known (ILev esque et al. 2010). Therefore, for the construc- 
tion of the afterglow SED, the GROND optical bands could be 
corrected for the contribution of the underlying host galaxy flux, 
but in the NIR only the maximum possible contribution of host 
galaxy light could be considered. 

Combining the GROND OBs 6 to 10 provides a good signal- 
to-noise ratio for the detection of the optical transient (afterglow 
plus underlying host galaxy) in g', r', i', z', J, and H (Table lA2l . 
Given the upper limits to the J, //-band magnitudes of the un- 
derlying host, the contribution of host galaxy flux in / and H 
at this time translates into an increase in brightness of the op- 
tical transient by at most 0.1 mag and 0.13 mag, respectively. 
To correct for this contribution, we added 0. 1 mag to the J as 
well as //-band data point and also increased the corresponding 
lcr error by 0. 1 mag. Performing the fit from g' to H then gives 
A) P t/NiR = 0.76 + 0.14 (Fig. 0, which is close to the observed 
mean for optical/NIR afterglows (about 0.6; cf. iGreiner et al.l 
l201lt iKann etai1[20Tol) . 

3.2. Multi-color light curve 



2 In the following we use the standard notation for the flux density 
of the afterglow, F v (t) oc r" v~ p . 

3 technical name for a pre-defined observing sequence 

4 http://www.sdss.org/dr7/algorithms/sdssUBVRITransform.html 



We combined TLS/GR ON D data with those oflXin et alJ (12010). 
ILevesque et al.l d2010l) . and lAntonelli et alJ d2009l) to obtain an 
/? c -band light curve composed of data sets published in refereed 
papers. Assuming a power-law SED of the afterglow and a non- 
evolving spectral index, the r'-band data were transformed into 
R c . For completeness, the V-band data from IXin et al. 
was also used and shifted into the R c band. The final /? f -band 
data set after the first bre ak at about 0.05 days ca n be fit using a 
single broken power law (Beuerma nn et al . 1999; Fig. |3). 

Fitting the data we find a late break in the light curve at 
around 0.4 days, in addition to a first break at around 0.02 days 
that had been previously known. This second break in the light 
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Fig. 3. /? f -band light curve of the GRB 090426 afterglow and 
best fit with a single broken power law after host galaxy sub- 
traction and correction for Galactic extinction. The fit starts 
at 0.05 days. For comparison, the g'i'z'JH bands (left y-axis; 
Tables lAll lA2l lA3l and the X-ray light curve (0.3 to 10 keV, 
Evan s et alj|2010t right y-axis) are also shown. 



curve was not seen in the previous data sets of GRB 090426. 
The g'r'i'z'JH GROND data show that this evolutionary phase 
is achromatic within the optical/NIR bands (Fig. [3). The most 
obvious interpretation is that this is a jet break. 

4. Discussion 

4. 1 . The light curve parameters 

Jet breaks are wi dely studied feat ures of long GRB afterglow 
light curv e s (cf . iFrail et aLll2001l) . In their long burst sample, 
IZeh et al.l (2006) find a nearly log-normal distribution of jet 
half-opening angles between 2 and 12 degrees, wit h the peak 
around 2 to 4 degrees (see also iRacusin et al.ll2.009h . For short 
GRBs, however, afterglow light curves ar e typically sparsel y 
sampled. Some cases seem uncollimated dGrupe et al.l 12006). 
while others display breaks and steep late slopes, which are evi- 
dence of c ollimation ( Burrow s et al ■l2006HMcBreen et al.l20ldt 
ISoderberg et al.ll2006l) . Most of these results, however, rely on 
the corresponding X-ray light curve. In the case of GRB 090426, 
on the other hand, there are basically no X-ray data available 
around the time of the late break in the optical light curve. The 
available data seem to indicate a smooth X-ray afterglow decay 
from 4000 s on up to the last X-ray detection at about 4 x 10 5 s 
dXin et al.ll2010l) . However, a break in the X-ray light curve at 
0.4 days, which is simultaneous with the break in the optical 
bands, is not ruled out. 

A satisfying fit of the entire optical and X-ray light 
curve can be obtained by assuming a two-component jet 
model dBerger et al.ll200l iPeng et alj2005t IRacusin et ai1l2008l: 
iFilgas et aLll201 it Fig. |4j). Within this framework, the observed 
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Fig. 4. Best fit of the optical afterglow using a two-component j e t 
model. Top: X-ray light curve (0.3 to 10 keV: lEvans et al.l20 10). 
Bottom: Combined /? e -band data set (the same as in Fig.[3j. 



afterglow light curve is the superposition of the radiation from 
two jets, a narrow and a wide jet. Even though it was not the 
aim of this paper to explore the validity of this model for GRB 
090426, we used it to fit the data and to shift X-ray and optical 
data points to the same time after the burst in order to obtain the 
SED from the optical to the X-ray band. 

According to the best fit, the narrow-jet component is de- 
scribed by a single broken power law with a\ - 0.48+0.04, = 
1.22 + 0.05, and a break time %\ = 290 + 20 s (while fixing the 
smoothness parameter n\ to 3). The second, wider component, 
follows a double broken power law with a\ = 0.46 + 0.15, as = 
2.43 ± 0.19 and break times t b2 = 9400 ± 3800 s (0.11 + 0.04 
days) and t b3 = 34500 ± 1800 s (0 39± .02 days; by requir- 
ing QT3 = -0.5 llPanaitescu & Kumail l2"bO01. «2=10, and «3=10; 
* 2 /d.o.f = 90.39/78 = 1.16). The optical and the X-ray light 
curve trace each other, implying that both are belonging to the 
same spectral regime. We caution, however, that after about 0.2 
days the X-ray data do not constrain the corresponding fit very 
much. 

The flat decay of the light c urve after its first break at t = fbi 
as well as the closure relations (Zhang & Meszaros 2004) show 
that this first break cannot be the jet break of the narrow-jet 
component. A plausible explanation of this part of the light 
curve could the c essation o f an energy in jection episode (cf. 
IZhangetalJ 12006: see also IXin et alJl20Toh . However, of pn- 
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mary interest here is only the late break of the optical light curve 
at t — ?b3. Regardless of the most suitable model, a single-jet 
model (Fig. [3} or a two-component jet model (Fig. 0), the de- 
duced break time as well as the post-break decay slope do not 
change. 



4.2. The closure relations and the jet parameter 

The Swift/XRT data dEvans et al.ll2010h show that for t > 10 4 s 
the spectral slope in the X-ray band is constant with fix - 
0.91 + 0.30. Combining this with the measured spectral slope 
in the optical/NIR bands (Sect. 13. Il l, and applying the closure 
relations, it follows that for t > 0.4 days a jet model is in reason- 
able agreement with the observations (Table lA~4l . Thereby, a jet 
model with si deways expansion is preferred. 

Following Sari et aT1dl999l) . the jet half-opening angle @ ISM 

jet 

for an ISM environment is ©] e s t M = ± (j^f' & {j^f^ ,where, 
£52 is the isotropic equivalent energy of the prompt emission 
in units of 10 52 erg, «o is the density of the ambient medium in 
cm -3 , rjy is the efficiency of the shock in converting the energy of 
the ejecta into gamma radiation, and % is the break time in days. 
The parameter a - is 1.0 for a single-jet scenario, while within the 
context of the two-component jet model the isotropic equivalent 
energy of the wide jet is around 10 % of the correspon ding num- 
ber for the narrow-jet compone nt dPeng et al. 1 120 05), i.e. it is 
k = 0.1. Using « = 10 cirT 3 dXin et al.ll2010l) and rj y = 0.2, 
for the observed brea k time at f/,= 0.4 days with £52 = 0.42^^ 
dLevesque et al.ll2010h . it follows that 0£ t = (6.5 + 0.4) degrees 

and E$ orr,w [1 to 10 4 keV] = (4.2 + 1.4) x 10 48 erg. Within the 
framework of the single-jet scenario, the break time is basically 
the same (Fig. [3j, k = 1, and it follows that jet = (4.8 + 0.3) 
degrees, as well as E c ™ = (2.3 + 0.8) x 10 49 erg. 

Before the suspected jet break at 0.4 days, the observed light 
curve decay between about 0. 1 and 0.4 days is quite shallow (0-4; 
Fig. |Ui, while the closure relations in this case, for a model with 
isotropic expansion, predict a steeper decay (a> 1). We caution, 
however, that the lack of data in this evolutionary phase makes 
it impossible to draw definite conclusions here. 



4.3. The afterglow compared to other short bursts 



Some authors have proposed (e.g., Zhan g et"ai1 12009|) a phe- 
nomenological classification based on the link with the progeni- 
tor, defining Type II in the case of a collapsar event and Type I in 
the case of merging compact objects, independent of the actual 

duration of the GRBs. 

Using the methods detailed in iKann et alJ ((2006), we cre- 
ated a composite light curve of the afterglow of G RB 09 0426 
and c ompared it to the afterglow samples of IKann et all 
(120101 Type II GRBs), and lKann et a"H d2008L Type I GRBs). 
Observationally, the afterglow is seen to lie in the faint end of 
the distribution of Type II GRB afterglows (Fig. IA.lt . especially 
in cases of late steep decay, and there are several Type I GRB 
afterg lows that are brig hter (e.g., GRB 0507 24, iBerger et al .1 
2005; GRB 05 1221 A ISoderberg et all [20061; GRB 060614 



Delia Valle et a 1l2006t iFvnbo et al.ll2006t iGal-Yam et al.l f2006: 
Man ganoetalJl2007h at this time. However, GRB 090426 lies 
at a much higher redshift (z = 2.6) than the aforementioned 
GRBs (z = 0.1 - 0.5). In the z = 1 system (Fig. E3, the af- 
terglow of GRB 090426 is clearly seen as an average afterglow 
relative to the Type II GRB afterglow sample. At t = 1 day 
(in the z = 1 frame), it would have had an absolute brightness 



of Mb = -22.02 ± 0.35, which is just one magnitude fainter 
(2cr) than the mean of the afterglow c omparison sample based 
on Swift detections dKann et alj 12010). On the other hand, the 
afterglow of GRB 090426 is at all times more luminous than 
any afterglow of the Type I sample exc ept for the controversial 
case of GRB 060121 dKann et alj|2008h . There is a strong indi- 
cation that, in spite of its very short duration, GRB 090426 is 
a Type II GRB, in accordance with other studies dZhang et al. 
2009; Levesaue et aDIIoTollXin et al.ll20Toh . 



5. Summary and conclusions 

We have presented TLS/GROND data of the optical/NIR af- 
terglow of GRB 090426, which show that the afterglow fea- 
tures a second break that was missed in all previously published 
data sets. On the basis of its achromaticity in the optical/NIR 
bands and in agreement with the closure relations, we have ar- 
gued that the late light curve break at 0.4 days is a jet break. 
Its calculated half-opening angle agrees well with the distribu- 
tion of half-opening angles found for long bursts. In addition, 
the observed luminosity of the afterglow also suggests that GRB 
090426 was related to a collapsar event. The interesting ques- 
tion then is whether the short duration of GRB 090426 in its 
host frame at z = 2.609 (Tgo = 0.33 s) can be explained within 
the framework of the collapsar model and how this compares 
to other long b ursts of similar short duration in their host frame 
dGreinerll201 ll) . More observational data of other short burst af- 
terglows are needed, not only to derive more reliable statistics 
but also to understand wether this short burst is an exception 
rather than the rule. 

In a recent paper jThdne et cd\\2011\) find further arguments that 
GRB 090426 was due to a collapsar event. 
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Appendix A: Observational Data and Afterglow 
Luminosity 



Table A.l. Log of the GROND observations (in case of the first 
epoch data these are OBs 1 to 10), with the magnitudes given 
in the AB system (not corrected f or Galactic ext i nction ). These 
results supercede the data given in lOlivares et al .1 (12009b . 
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21.46 (08) 


21.30 (08) 
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>19.8 




46268 


22.03 (04) 


21.67 (05) 


21.59 (09) 


21.41 (11) 


>21.6 


>20.8 


>19.8 




47037 


22.00 (05) 


21.82 (05) 


21.55 (08) 


21.47 (14) 


>21.6 


> 21.0 


>19.9 




47812 


22.08 (05) 


21.83 (04) 


21.64 (08) 


21.29(15) 


>21.7 


> 21.1 


>20.2 




49112 


22.17 (03) 


21.84 (03) 


21.67 (05) 


21.50 (07) 


>22.1 


> 21.4 


>20.4 




50930 


22.28 (03) 


21.96 (03) 


21.79 (06) 


21.52(07) 


>22.1 


> 21.3 


>20.4 




52756 


22.34 (03) 


22.10(04) 


21.87 (05) 


21.71 (09) 


>l£,l 


> 21. 3 


>20.J 




54571 


22.44 (04) 


22.20 (03) 


21.96 (07) 


21.86 (05) 


>22.0 


> 21.2 


>20.5 




56374 


22.58 (04) 


22.29 (04) 


22.07 (06) 


21.93 (11) 


>22.2 


> 21.6 


>20.4 




139787 


24.03 (10) 


23.80(11) 


23.63 (17) 


23.38 (18) 


>22.7 


>21.9 


>21.1 




222822 


24.23 (11) 


24.24(12) 


23.84 (15) 


>24.1 


>22.7 


>21.8 


>21.3 




Table A.2. 


Log of the GROND observations for the combined 










OBs 1 to 5 and OBs 6 to 10, with the magnitudes ; 


given in the 










AB system. 


Data are not corrected for Galactic extinction. 










Time (s) 


g 1 


r> 






J 


H 




K 


46268 


22.05 (03) 


21.74 (03) 


21.65 (05) 


21.44 (07) 


20.96 (02) 


20.92 (20) 


> 20.1 


52723 


22.28 (02) 


21.99 (02) 


21.77 (03) 


21.62 (04) 


21.39 (09) 


21.23 (12) 


>20.6 



Table A.3. Log of the TLS observations, given in the Vega sys- 
tem. Data are not corrected for Galactic extinction. 



Time(s) 




h 


26868 


20.90 (13) 




27146 


20.91 (14) 




27607 


21.05(19) 




28063 


20.97 (14) 




26263 




20.48 (18) 
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Table A.4. Predicted temporal decay slopes a for t > 0.4 days 
for various afterglow scenarios based on the measured spectral 
slopes /J opt = 0.76 ± 0.14 and f3 x = 0.91 + 0.30. 



Model Optical 


X-rays 


■S/o"s,opI 




Isotropic case 










ISM, v < v c 1.13 ±0.21 


1.37 ± 


0.45 


-4.61 


-2.18 


ISM, v > v c 0.63 ± 0.21 


0.87 ± 


0.45 


-6.37 


-3.20 


Wind, v<v c 1.63 ±0.21 


1.87 ± 


0.45 


-2.84 


-1.16 


Wind, v>v c 0.63 ± 0.21 


0.87 ± 


0.45 


-6.37 


-3.20 


Jet with sideways expansion 


v < v c 2.50 ± 0.28 


2.82 ± 


0.60 


0.21 


0.62 


v>v c 1.50 ±0.28 


1.82 ± 


0.60 


-2.75 


-0.97 


Jet without sideways expansion 








v<v c 1.88 ±0.21 


2.12 ± 


0.45 


-1.96 


-0.64 


v>v c 1.38 ±0.21 


1.62 ± 


0.45 


-3.73 


-1.67 



Columns 4 and 5 give the difference between the predicted and the 
observed (a = 2.43 ± 0.19; the parameter a s in Sect. I4.lt temporal 
decay slope, normalized to the square root of the sum of their quadratic 
errors, with s = (a pledicted - (Wved), cr] = o- 1 ^^ + <T^ bs(:rved . 
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1E-4 1E-3 0.01 0.1 1 10 100 

t (days after burst in the observer frame) 



Fig. A.l. The afterglow of GRB 090426 (thick blue line) in comparison with the afterglows of Type II (thing gray lines) and Type I 
(red symbols and lines; squares c onnected by spline s are d etections, downward triangles connected by thin dashed lines are upper 
limits) GRBs from the sample of iKann et all (120081 l2010h . These afterglows have been corrected for Galactic extinction, and the 
host galaxy contribution has been subtracted where possible (also in the case of GRB 090426). The afterglow of GRB 090426 is 
seen to be among the faint Type II GRB afterglows, but it is brighter than most Type I GRB afterglows or limits thereon. 
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t (days after burst in the observer frame assuming z = 1) 



Fig. A.2. The afterglow of GRB 090426 after it has been shifted to the z — 1 system, again in comparison with the samples (shifted 
using the same method) of iKann et alj d2008ll20lol) . The labelling is identical to that in Fig. IA.ll It can now clearly be seen that the 
afterglow of GRB 090426 is readily comparable to the afterglow of Type II GRBs (collapsar events), while it is much brighter than 
any Type I GRB afterglow (merging compact objects), with the exception of GRB 060121, which is a controversial case. 
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Fig. A.3. X-ray to optical spectral energy distribution of the afterglow of GRB 090426 at t — 8ks after the burst, calculated based on 
the fitted light curve (Fig. |4j». The fit uses A^ al = 0.015 x 10 22 cirT 2 and corresponds to a negligible host extinction, a gas column 
density of A^ ost = 0.46°-," x 10 22 cirT 2 , and a spectral slope of /fox = 0.90 + 0.03 (x 1 = 10.85 with 7 d.o.f.,^ v = 1.55). 



